Strong enhancement of extremely energetic proton production in central heavy ion 

collisions at intermediate energy 
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The energetic proton emission has been investigated as a 
function of the reaction centrality for the system ^*Ni + ^*Ni 
at 30A MeV. Extremely energetic protons {E^^ > ISOMeV) 
were measured and their multiplicity is found to increase al- 
most quadratically with the number of participant nucleons 
thus indicating the onset of a mechanism beyond one and 
two-body dynamics. 

PACS numbers: 25.70.-z, 25.75.Dw 



Heavy ion collisions at intermediate energy allow to 
investigate the properties of nuclear matter far from sta- 
bility. Dynamical calculations show that, in the early non 
equilibrated stage of the reaction, high temperatures and 
densities are reached. Since heavy ion reactions at inter- 
mediate energy are described in terms of mean field and 
two body collisions, experimental observables sensitive 
to basic ingredients of the models such as the nucleon- 
nucleon (NN) cross section in matter and the mean field 
potential are needed to probe the nuclear dynamics and 
the equation of state of nuclear matter (EOS). In partic- 
ular, particles such as subthreshold mesons or energetic 
photons and nucleons are expected to provide informa- 
tion on the nuclear dynamics at the pre-equilibrium stage 
10,^ and reference therin. Experimentally the presence 
of a pre-equilibrium component in the light particle emis- 
sion has been observed in the energy spectra The 
impact parameter dependence of the pre-equilibrium pro- 
tons, investigated for several systems, shows that the av- 
erage energetic proton multiplicity increases almost lin- 
early with the number of participant nucleons . These 
features support the hypothesis of pre-equilibrium proton 
emission due to incoherent quasi- free NN collisions 
Moreover, experimental evidences, such as the observed 
7-proton anti-correlation and the energetic proton an- 
gular distributions reported in reference provide in- 
formation on the space-time origin of the energetic pro- 



tons indicating that a relevant fraction is emitted from 
first chance NN collisions in the interaction zone. These 
evidences show that the energetic protons are emitted in 
the first stage of the reaction according to expectations 
and therefore are good candidates to probe the pre- 
equilibrium phase. On the other hand, the observation of 
extremely energetic nucleons or deep subthreshold parti- 
cles over a broad range of incident energy addresses the 
question of which mechanisms could enable to concen- 
trate a relevant fraction of the available energy in the 
production of a single energetic or massive particle [^j. 
In fact, the emission of particles with energy or mass 
much larger than that provided by the coupling of the 
nucleon Fermi motion with the relative motion of the 
colliding nuclei is a challenging aspect of heavy ion col- 
lisions both experimentally and theoretically, due to the 
very low production rates and the lack of information on 
the production mechanism. 

In this letter we present results concerning the emis- 
sion of protons with energy extending up to almost 20% 
of the total available energy in the reaction ^^Ni + ^^Ni at 
30A MeV. Since these energies largely exceed the maxi- 
mum energy expected in first chance NN collisions due to 
the coupling of the relative motion with a sharp nucleon 
Fermi momentum distribution (kinematical limit), this 
investigation can also provide a clue for the comprehen- 
sion of the deep subthreshold particle emission in this en- 
ergy domain. In particular, for the first time, the proton 
multiplicity as a function of the impact parameter was 
measured for extremely energetic protons {E^^ > 130 
MeV). A strong non linear dependence on the number 
of participant nucleons is observed thus providing im- 
portant information on the production mechanism. A 
detailed comparison with a microscopic transport model 
was also performed aiming to extract information on the 
dynamics at pre-equilibrium and the nature of energetic 
proton emission. 
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The experiment was performed at Laboratori Nazionali 
del Sud with the MEDEA § and MULTICS (lO| appara- 
tus. A ^^Ni beam at 30A MeV dehvered by the Tandem 
and Superconductive Cyclotron (CS) acceleration system 
bombarded a ^^Ni target 2 mg/cm^ thick. MEDEA con- 
sists of a ball made of 180 BaF2 detectors placed at 22 cm 
from the target which covers the polar angles from 30° 
to 170°. The BaF2 permits to detect and identify LCP 
(Ep < 300 MeV) and photons up to E-^ « 200 MeV. The 
time of flight and pulse shape discrimination analysis al- 
lows to clearly identify photons, protons, deuterons and 
tritons, alpha particles [^. The response of BaF2 crys- 
tals to LCP has been investigated using monoenergetic 
particle beams and a calibration procedure based on the 
7 calibration with 7 sources and cosmic rays has been 
established as described in ref. Q . The MULTICS ar- 
ray is made of 55 telescopes covering the angular range 
3° < 9iab < 28°. Each telescope consists of an Ionization 
Chamber, a Silicon detector and a Csl crystal, and allows 
the identification of charged particles up to Z = 83. The 
threshold for charge identification was about 1.5 A MeV 
pof . The total geometric acceptance was greater than 
90% of An. 

Energetic protons were detected in coincidence with 
photons, light charged particles {Z = 1,2) (LCP) and 
intermediate and heavy fragments on an event by event 
basis, thus allowing a rather complete description of the 
reaction dynamics as well as an estimate of impact pa- 
rameter which represents a crucial problem in this energy 
domain. Due to the very low cross section expected for 
extremely energetic protons, high statistics spectra are 
needed. To increase the fraction of events containing en- 
ergetic protons, the main trigger required the presence 
of at least one BaF2 signal above a threshold level cor- 
responding to proton energy of about 30 MeV. Moreover 
the coincidence with MULTICS reduced the cosmic ray 
contamination to a negligible level. All the MEDEA de- 
tector with 6 > 75° and a few detectors of the forward 
rings took part in this trigger. Events corresponding 
to a minimum bias trigger, defined by the OR between 
MEDEA and MULTICS, were also scaled down and reg- 
istered. Altogether approximately 4 x 10* events have 
been collected and analysed. 

According to the standard three moving source analy- 
sis, the high energy proton emission at large polar angles 
can be described by a source emitting with velocity close 
to the half beam velocity and a high inverse slope param- 
eter. A selection of the energetic protons emitted from 
this intermediate velocity source is possible by applying 
kinematical constraints (i^^'"'' > 40 MeV, 6'iab > 42°) 
| p2[ . The experimental proton spectra, transformed in 
the NN frame {v/c — 0.127) are reported in fig. la (full 
symbols) together with the intermediate source compo- 
nent of the fit (solid lines) for inclusive data. The inverse 
slope parameter deduced from the maxwellian fit with 
a volume emission is in good agreement with the sys- 



tematics (T « 11 MeV) Very energetic protons are 
observed in the spectra, with energy well above the kine- 
matical limit expected in the hypothesis of first chance 
NN collisions and sharp Fermi momentum distribution 

cam 

) (arrows of fig. 1). 
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FIG. 1. (a) Experimental inclusive proton multiplicity 
spectra in the NN reference frame at different polar angles 
(87°, 110°, 113°, 126°, 140°). (b) Experimental and BNV pro- 
ton energy multiplicity spectrum in the NN reference frame 
for central collisions (< b >^ 2.5 fm, 98° < 6'ivjv < 124°). Ex- 
perimental data (full squares), local Skyrme interaction (open 
squares) and momentum dependent interaction (open circles) 
are reported. The arrows indicate the expected kinematical 
limit for first chance NN collisions. 

Extremely energetic protons were already observed in 
heavy ion collisions (see and refs. therein), but those 
experiments did not allow a conclusive answer on the pro- 
duction mechanism to be drawn. In order to explain the 
production of particles with energy or mass larger than 
the energy available in NN collisions including the Fermi 
motion, several hypotheses, such as the presence of high 
momentum tails in the nucleon momentum distribution, 
the onset of multistep processes, cooperative effects, the 
presence of dynamical fluctuations, have been proposed. 

To gather more information on the energetic proton 
emission, we have performed simulations with a BNV 
code which is based on mean field and two body dynamics 
p3[ . To explore the sensitivity to basic ingredients of the 
calculations two different potentials were implemented in 
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the BNV code: one using a local Skyrme interaction for 
the mean field (open squares of fig. lb) and another using 
a Gale-Bertsh-Das Gupta momentum dependent inter- 
action (GBD) (open circles of fig. lb) ||l^. Concerning 
the reaction dynamics, both calculations exhibit similar 
features leading to the formation of a heavy residue in 
central collisions, whereas binary collisions dominate at 
larger impact parameters. However, although momen- 
tum dependent effects are expected to be more impor- 
tant at higher bombarding energy, it is interesting to un- 
derstand how far their inclusion can affect the energetic 
proton production which is ruled by a delicate balance 
between the mean field and the nucleon-nucleon cross 
section. The experimental and calculated spectra in the 
NN frame are shown in fig. lb for central collisions. To 
achieve a satisfactory statistics for energetic proton pro- 
duction, about 250 events have been simulated per im- 
pact parameter (using 200 test particles per nucleon). 
The criteria for impact parameter selection and assign- 
ment will be discussed in detail in the following. Our 
results point out remarkable differences: the spectrum 
calculated with the local Skyrme interaction strongly un- 
dershoots the data and exhibits a lower inverse slope pa- 
rameter than the experimental one, while the spectrum 
calculated with the momentum dependence is in better 
agreement with the data concerning both the yield and 
the slope (at least up to « 110 MeV). Due to the less at- 
tractive mean field, we observe in the GBD calculations a 
larger fraction of escaping particles. For the same reason, 
these particles can also be more energetic. Moreover, the 
calculated yields should be slightly reduced since in this 
kind of calculation only free nucleons are emitted while 
in the reaction also complex particles are emitted and 
observed experimentally 

From the comparison with the experimental data one 
can get a deeper insight into the behaviour of nuclear 
matter at large density and temperature. Therefore, with 
the aim of improving the overall understanding of the 
energetic proton emission and disentangling between the 
various hypotheses for the production of the most ener- 
getic protons, we have investigated the impact parameter 
dependence. Indeed, the dependence of multiplicity on 
the number of nucleons participating in the reaction can 
provide information about a change in the production 
mechanism. A stronger than linear increase of the multi- 
plicity as a function of the number of participant nucleons 
has been observed, at much higher incident energy, in the 
deep subthreshold production of K+ |l^ , rj [|l6| and high 
transverse- mass tt*^ ||^. 

At energy as low as 30A MeV, the fluctuations on 
global variables such as the charge particle multiplicity 
and transverse energy affect the determination of the im- 
pact parameter especially for the most central collisions 
. To solve the problem of the event selection over a 
wide range of impact parameters, we exploit the reaction 
mechanism and hard photon multiplicity information to 



determine the size of the interaction zone ||l^,|T^. In- 
deed, the detection of heavy fragments from projectile- 
like fragments to evaporation residues and their relative 
velocities allows to select classes of events with differ- 
ent centrality. In particular, the most central collisions 
were selected requiring the presence of an evaporation 
residue with velocity close to the centre of mass veloc- 
ity and charge higher than the projectile charge. On the 
other hand, the detection of fragments originating from 
projectile fragmentation and deep inelastic collisions was 
exploited to reach classes of events spanning the range of 
impact parameters from peripheral to mid-central. The 
same procedure was also applied to some classes of events 
selected in terms of total charged particle multiplicity. 
Finally, to obtain a quantitative estimate of the impact 
parameter, the hard photon multiplicity has been cal- 
culated for the various classes of events. The number 
of participant nucleons Apart {b) has been extracted from 
the hard photon {E^ > 30MeV) multiplicity, according 
to the relation M^{b) = ■ Nnp{h) ~ • 0.5 • Apart{b) 
where P-^ is the probability of emitting a hard photon 
in a np collision deduced from inclusive data {P^{Ej > 
aOMey) w 2.7 ■ 10~5) and 7V„p(6) is the number of first 
chance np collisions occurring in the overlap region. This 
relation has been satisfied by several experiments which 
show that the hard photon multiplicity provides a snap- 
shot of the participant region ||20|-p^ . 

In fig. 2 the average proton multiplicity is reported as a 
function of the number of nucleons participating Apart (b) 
in the collision for different energy bins in the NN ref- 
erence frame (60 ^ 80 MeV (Mp(60)), 100 ^ 120 MeV 
(Afp(lOO), 130-150 MeV (Mp(130)). The experimental 
proton multiplicity (full squares) displays the expected 
linear dependence on Apartib) for energy close to the 
kinematical hmit (60 < Ep < 80 MeV fig. 2a), while a 
stronger dependence is observed with increasing proton 
energy, in particular the multiplicity of extremely ener- 
getic protons (Mp(130)) exhibits an almost quadratic in- 
crease with Apart (fig- 2c). 

The BNV calculations, filtered with the experimental 
apparatus, are also reported in fig. 2 (open circles). The 
Apart{b) assignment relies on the hypothesis of a geo- 
metrical correlation between b and Apart{b) The 
calculations have been scaled by a factor 0.6 to allow a 
better comparison with the data. This scaling is con- 
sistent with the yield reduction needed to account for 
complex particle emission. Within this assumption, a 
good agreement with the data is observed in fig. 2a and 
fig. 2b, confirming that the energetic proton production 
is described with good accuracy up to « 110 MeV. On 
the other hand, BNV calculations fail in the most ener- 
getic bin (Mp(130), fig. 2c) where the almost quadratic 
dependence on Apart observed experimentally is not re- 
produced thus showing the onset of effects beyond the 
mean field and two body collisions. It is interesting to 
notice that a non linear dependence is observed, both 
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experimentally and theoretically, also for Mp(lOO). The 
calculation can account for this behaviour due to the in- 
creasing importance of multistep two-body collisions in 
the production mechanism of protons with energy higher 
than the kinematical limit |^. However, this mechanism 
seems not to be able to explain the almost quadratic be- 
haviour observed for Mp(130). Indeed, for extremely en- 
ergetic protons, this multistep process is associated with 
larger time scale. Therefore the system can emit nucle- 
ons and rapidly evolves far from the initial geometrical 
overlap configuration. This can explain the weaker de- 
pendence on the impact parameter observed in the cal- 
culations (fig. 2c). 
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FIG. 2. Proton multiplicities as a function of the number 
of participant nucleons (see text) for different energy bins. 
Experimental values (solid squares) and momentum depen- 
dent BNV calculations (open circles) scaled by a factor 0.6 
are reported. Only protons emitted in the angular range 
75° < 9iab < 138° are considered. To compare the trends 
a linear and a quadratic dependence are reported (solid lines) 
in panels a) and c), respectively. 

The observed behaviour of the multiplicity of very en- 
ergetic protons on the number of participant nucleons 



(fig. 2c) puts constraints on the mechanism responsi- 
ble for the production of extremely energetic protons. 
Dynamical fluctuations |^ are not expected to lead to 
the Apart quadratic behaviour observed experimentally. 
Other effects, such as high momentum tails, are weakly 
dependent on density and, at the energy considered, den- 
sity variation from central to peripheral impact parame- 
ters are small [Q. Cooperative effects, where more nu- 
cleons or clusters of nucleons participate in the collision, 
seem very promising and should be investigated. 

In summary, the energetic proton production has been 
investigated up to proton energy corresponding to about 
20% of the total energy available in the system. The 
energetic protons up to « 110 MeV are emitted as a 
consequence of NN collisions in the first stage of the re- 
action and their characteristics are well reproduced by 
BNV calculations which include the momentum depen- 
dence in the effective potential. On the other hand, the 
BNV approach fails to explain the almost quadratic de- 
pendence on the number of participant nucleons of the 
yield of very energetic protons {E^^ > 130MeV). This 
behaviour calls for the introduction of mechanisms be- 
yond the mean field and two body nucleon-nucleon colli- 
sions such as cooperative effects. These results shed some 
light on the emission of extremely energetic protons and 
can improve the understanding of the mechanism respon- 
sible for deep subthreshold particle production. 

We thank Prof. M. Di Toro for the critical reading 
of the manuscript and the LNS staff for the high-quality 
beam and the support during the experiment. 
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